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ON  COMPUTING  THE  NORM  OF  A  TRANSFER 

MATRIX 


S.  Boyd,  V.  Balakrishnan  * 
Information  Systems  Laboratory 
Electrical  Engineering  Department 
Stanford  University 

Abstract 

We  present  a  simple  bisection  algorithm  to  com¬ 
pute  the  Hqo  norm  of  a  transfer  matrix.  The 
bisection  method  is  far  more  efficient  than  algo¬ 
rithms  which  involve  a  search  over  frequencies, 
and  moreover  can  compute  the  Hoo  norm  with 
guaranteed  accuracy. 

1  Preliminaries 

Throughout  this  paper  A,  B,  C,  D  will  be  real 
matrices  of  sizes  n  x  n,  n  X  m,  p  x  n,  and  p  X  m, 
respectively.  We  refer  to  the  linear  dynamical 
system 

x  =  Ax  +  Bu  ,  . 

y  =Cx  +  Du  ’ 

as  the  system  {A,B,C,D}.  We  refer  to  H(s)  — 
C{sl  —  A)-1  B  +  D  as  the  transfer  matrix  of  the 
system  {A,B,C,D}. 

A  is  stable  means  that  all  eigenvalues  of  A  have 
negative  real  part.  If  A  is  stable,  we  define  the 
H^,  norm  of  the  transfer  matrix  H(s)  to  be 

Halloo  =  8UPRe»>o£r"»«*(^(,s))  (2) 

=  SUPu*R  K  * 

where  <7mo*(-)  denotes  the  maximum  singular 
value  of  a  matrix,  that  is,  crmax(F)  =  AJ^(F*F). 

The  Hqo  norm  of  a  transfer  matrix  arises  often 
in  control  theory.  An  important  interpretation 
of  j|/f  ||oo  >s  88  the  Lj  or  RMS  gain  of  the  system 
(1)  (see  e.g.  [2]). 
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]|//||oo  is  usually  ‘computed’  by  searching  for 
the  maximum  of  ))  over  u;  €  R.  Ob¬ 

vious  problems  associated  with  such  a  method 
are  (a)  determining  the  range  and  spacing  of  the 
frequencies  to  be  checked,  and  (b)  the  large  num¬ 
ber  of  computations  involved  (a  si ngular\ value 
decomposition  (SVD)  is  often  performed  at  each 
frequency  point).  The  problem  (a)  is  particularly 
evident  when  A  has  lightly  damped  eigenvalues. 

We  propose  instead  a  bisection  method  in¬ 
spired  by  Byers’  bisection  method  for  measuring 
the  distance  of  a  stable  matrix  to  the  unstable 
matrices  [3).  The  bisection  method  not  only  in¬ 
volves  less  computation,  but  has  the  advantage 
of  computing  |j//||oo  with  a  guaranteed  accuracy. 

2  SVs  of  H  via  a  Hamilto¬ 
nian 

We  start  by  establishing  a  connection  between 
the  singular  values  of  the  transfer  matrix  and  the 
imaginary  eigenvalues  of  a  certain  Hamiltonian 
matrix.  Let  7  >  0,  and  not  a  singular  value  of 
D.  Define 

M7  = 

A-BR~1DtC  -jBR-'B1, 

7CrS-1C  -At  +  CtDR-'Bt 

where  R  =  (DTD  -  72/)  and  S  =  (DDT  -  72/). 
Note  that  M7  is  a  Hamiltonian  matrix. 

The  following  theorem  relates  the  singular  val¬ 
ues  of  H(ju)  and  the  imaginary  eigenvalues  of 
A/7. 
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Theorem  X  Assume  A  has  no  imaginary  eigen¬ 
values,  7  >  0  is  not  a  singular  value  of  D,  and 
u-'o  G  IX. 

Then,  7  is  a  singular  value  of  II(jaj0)  <=> 
(My  —  jio0I)  is  singular. 

Remark  1: 

There  ar-  no  observability,  controllability,  or 
stability  conditions  on  the  system  {A,  B,C,  D}. 
A  simple  consequence  of  Theorem  1  is 

Theorem  2  Let  A  be  stable  and  7  >  <rmax(D). 
Then  ||//i|«,  >  7  4=>  My  has  imaginary  eigen¬ 
values  (i.e.  at  least  one). 

Remark  2: 

The  imaginary  eigenvalues  of  M^j^  are  ex¬ 
actly  the  frequencies  for  which  crmax(II(j u>))  = 

Remark  3: 

Theorem  2  is  also  readily  derived  via  several 
methods,  e.g.,  quadratic  optimal  control  [4]  or 
spectral  factorization  [5]. 

3  Bisection  Algorithm 

Theorem  2  suggests  a  bisection  algorithm  for 
computing  |j//||oo-  Let  jn  and  7^  be  some  lower 
and  upper  bounds,  respectively,  on  ||//||oo-  For 
example,  one  could  use  the  bounds  derived  by 
Enns  and  Glover, 

71*  =  max{crmar(JD),  cr/n} 
n 

7ui  =  &max  (D)  +  2  (Tjji 

txl 

where  erm  are  the  Hankel  singular  values  of  the 
system  {A,  B,  C ,  D }  [6,  7]. 

The  bisection  algorithm  is  as  follows: 

U  ■=  7/41 
7 U  '■=  7u4i 
repeat  { 

7  :=  (7 L  +  7w)/2; 

Form  Af7; 

if  My  has  no  imag.  eigenvalues,  7#  :=  7, 
else  7*.  =  7;  } 
until  {  7ff  ~  7l  <  2e  yL). 

Note  that  we  always  have  7 1  <  ||^f||oo  <  7 //. 
On  exit,  (7^  +  7h)/2  is  guaranteed  to  approxi¬ 
mate  ||//||oo  within  a  relative  accuracy  of  e,  i.e., 

l(7L  +  7ff)/2-||i/||oo|<e||//|U 


Remark  4 [: 

Checking  if  My  has  any  imaginary  eigenvalues 
can  be  done  in  a  finite  number  of  steps  via  a 
Sturm  sequence  test  on  the  characteristic  poly¬ 
nomial  [8]. 
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